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Abstract. This article addresses the problem of enumerating the tilings 
of a plane by lozenges, under the restriction that these tilings be doubly 
periodic. Kasteleyn's Pfaffian method is applied to compute the gener- 
ating function of those permutations. The monomials of this function 
represent the different types of tilings, grouping them according to the 
number of lozenges in each orientation. We present an alternative ap- 
proach to compute these types. Finally, two additional classes of tilings 
are proposed as open enumeration problems. 

1 Introduction 

We consider the tiling of the plane by equilateral triangles, assuming that the 
vertices of those are the points of the lattice A spanned by u = (1,0), v = 
(1/2,-^/3/2). 




Figure 1: Equilateral triangle tiling 



By merging two adjacent triangles of this tiling we obtain a rhombus or 
lozenge. It is obvious then that the plane can also be tiled with such lozenges. 
These tiles may take three orientations, which we denote by the notation set in 
Figure [2j 




Figure 2 



There is a rich literature devoted to studying the different lozenge arrange- 
ments that tile certain bounded regions of the plane (see [10, 2, 4, 9 and the ref- 
erences therein) . The simplest scenario appears in the case of a convex hexagon 



Figure 3: A lozenge tiling in a semiregular hexagon 



drawn over the grid of Figure [T] with pairs of opposite sides of the same length 
(a semiregular hexagon). Its lozenge tilings are associated by an appealing re- 
jection with the plane partitions whose parts are bounded by the lengths of the 
hexagon sides. We recall that the solid Young diagram of plane partition is an 
arrangement of unit cubes located in the positive octant of K 3 and satisfying 
that if there is a cube at position k), a < i, b < j, and c < k, then position 
(a, b, c) is neither empty. 

A classical formula (proved in llj) states that the number of plane partitions 
fitting in a box of sides a, b, and c is: 

a—1 b—1 c-1 . , _ 

Throughout this article, we employ bold letters to denote vectors a = 
(01,02), considered as columns when using matrix notation: a = (01,02)*. We 
are interested in lozenge tilings of the whole plane which are doubly periodic. 
More explicitly, let 

?2x2 



B = 



at 


b 2 ' 




62 . 



be a 2-rank matrix and consider the sublattice 

A = [u|v]BZ 2 C A . (1) 

The tilings wc arc interested of arc invariable under the translations by aiU+a2V 
and 61U + &2V; or equivalently, by any element of A. Dealing with infinite tilings, 
unlike in the finite regions case, we do not have a boundary to "support" our 
arguments. 

We will define the type of a tiling as the number of lozenges in each orien- 
tation it consists of. Note that, as can be checked using the argument exposed 
in [3J; in the mentioned case of a semiregular hexagon with sides a, 0, and c, all 
its lozenge tilings have a constant number of lozenges in each orientation (afe, 
oc, and be). 

As we have said, a lozenge tiling is built merging pairs of triangles under 
certain restrictions. Therefore, lozenge tilings correspond to perfect matchings 
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Figure 4: A doubly periodic tiling 



in an associated graph (see Figure [8]). Perfect matchings enumeration is a 
rich field of research (see [13 for an illustrating survey), stimulated by several 
problems in the domain of Physics and Chemistry. 

Kasteleyn developed a method for enumerating the perfect matchings of a 
planar graph by means of a Pfafhan computation [6 , 7 , 8J . In the first of these 
references, the problem is also solved for a non-planar graph, which can be 
embedded in a torus. This is accomplished by computing a linear combination 
of four Pfafhans. In Section [3J we explain how this methods applies to counting 
doubly periodic lozenge tilings. The output is a polynomial generating function 
Z(L, D, R) 6 Z[L, D, R], where the coefficient in UD 3 R k equals the number of 
tilings of type k). 

In Section [4j we introduce a different approach to compute the different 
types of doubly periodic tilings modulo a given period A. We prove that the 
pair of height increments in the (infinite analogue to the) solid Young diagram 
through a pair of vectors spanning A characterises the type of the tiling. As a 
consequence, the possible types correspond to lattice points in a certain triangle 
associated to a lattice basis B, in the sense of 0. The vertices of this triangle 
correspond to the three uniform tilings (those with all of its lozenges arranged 
in a constant orientation). 

Finally, in Section [5j we propose two ways of grouping together significantly 
similar tilings, whose enumeration remains, up to our knowledge, an open prob- 
lem. Before all, let us fix the notation and definitions we will need: 

2 Definitions 

There are two types of triangles in Figure [Tj upwards and downwards-pointing. 
Any lozenge contains a triangle of each type. In this article, we refer to generic 
elements of these classes by the symbols A and V . We identify both sets of 
triangles with Ao by means of the following convention: a point in Ao represents 
the right-most upwards-pointing and downwards-pointing triangles which have 
that point as vertex. 

Let A be defined by Equation |l]). We refer to the index of A in A simply 
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as the index of A: 

[A : A] = |det£| = 2(volA)/V3. 

In order to define a A-periodic tiling, we need to decide with which of the three 
adjacent V is merged every A in a fundamental set whose size is the index 
of A. The following map returns the difference within a lozenge on its 

orientation as input. 

£: {L, D, R} ->■ 
L i y 
D i-> 
R i-> 

Definition 1 Let A C A be a 2-rank lattice. We define a A-periodic tiling 

as a map 

r:A ^{L,D,R} 
satisfying the following two axioms: 
i) Compatibility: 

Vx e A Vy e A, r(x) = r(x + y). 

ii) Tiling: The following map is bijective: 

f : A -> A 

x i } x + £(r(x)). 

We use the notation T\ for the set of A-periodic tilings. 
Note that the second axiom is equivalent to (see Figure [5]): 
ii') Vx £ Aq, exactly one of the following conditions is satisfied: 
t(x) = R, t(x + u)= L, r(x + v) = D. 




Figure 5: Axiom ii) Exactly one of the three possible lozenges occurs. 

As we have already mentioned, it is enough to define r on a representative of 
each class modulo A. We consider the induced mapping: 

f : A /A ->■ {L, D, R} 
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and define the type of a A-periodic tiling as: 

i(r,A)= (#f- 1 (L),#f- 1 (£>),#f- 1 (fl)) £N 3 

Infinite lozenge tilings can be represented by an infinite analogue to a solid 
Young diagram, namely, the complement of a subset of Z 3 closed under addition 
of elements in the semigroup N 3 . In other words, a staircase diagram in three 
dimensions; with two possible identifications. As Figure k5] shows, there are two 
lozenge tilings of the unit hexagon. One of them is prominent (a solid cube), 
and the other one is a "hole" , limited by three walls. In the rest of this article, 
we consider that the left design of Figure [6] is the solid cube. 



Figure 6: An outside and an inside corner. 



With this identification, the vertices of a lozenge tiling can be labelled by 
three-dimensional coordinates and given a height function, defined by the sum 
of these coordinates. The application of height labels in this context dates back 
to 1, 15 and has been useful in the study of tilings of bounded regions. 

An edge in the tiling is associated with a coordinates and a height increment: 





F 


/ 


u 


(-1,0,0) 


-1 


— u 


(+1,0,0) 


+1 


V 


(0,0, +1) 


+1 


— V 


(0,0,-1) 


-1 


u — V 


(0,+l,0) 


+1 


v — u 


(0,-1,0) 


-1 




Note that this effect of an edge on the height of a vertex is only valid is 
the edge does appear in the lozenge tiling. For instance, if two vertices x and 
x + u are not connected by u, then r(x) = D and they are connected by the 
concatenation (v, u — v). Therefore, the coordinates increment is not (—1,0,0), 
but (0, +1, +1). In general, if the height increment of an existing edge is h, the 
effect of the same non-occurring edge is — 2h. 

Definition 2 Let A C Ao be a 2-rank lattice and r an A-periodic tiling. We 
define a path in r as a succession (Xi)fL C Ao satisfying 

Xj+i - Xj e {±u, ±v, ±(u - v)} 



and: 
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• x i+ i Xj = u =>• r(xj) 7^ L> 

• x i+ i - Xj = v r(xj) 7^ L 

• x i+ i - Xj = v - u =4> r(x 4 - u) 7^ _R 

• x i+ i Xj = u r(xj - u) 7^ D 

• x i+1 - x 4 = -v t(x 4 - v) L 

• Xj + i Xi = u v =>■ r(xi — v) 7^ i? 

Each path p = (xi)^ in t can be associated with a height increment, in the 
following way: 

N-l 

h r(p) = /( X *+l ~ X *)' 

i=0 

For a given tiling t, the height increment of a path only depends on its 
extreme points. Therefore, setting h T (0) = 0, we can associate a height to every 
point in Ao, defining a mapping h T over Ao. As we usually treat elements of 
A by its coordinates on basis (u|v), we will employ the following change of 
coordinates: 

e T : Z 2 — > Z 

a i y h T ({u\v]a). 




Figure 7: The h function. 



Note that as a lozenge tiling is the projection of an pile of cubes along lines 
parallel to vector (1, 1, 1), the height h T (x) of a point (together with its projected 
position in the plane) determines its three-dimensional coordinates, which by 
analogy we may denote i? r (x). 



G 



Figure 8: The bipartite graph associated to a sublattice of Aq 



Proposition 3 Let A C Aq be a 2-rank lattice, t a A-periodic lattice and x = 
A1U + A2VGA0. Then, 

H T (x) = (-A 1)0 , A 2 ) + Mx)-A, + A 1 (1| 1; iy 

Proof. The maybe improper "path" (0, u, 2u, . . . , Aiu, Aiu + v, . . . x) would 
give a height label h = A2 — Ai and coordinates label H = (— Ai,0, A2). Any 
wrong step can be replaced by a concatenation of two, changing the height by 
±3 and the coordinates by ±(1, 1. 1). ■ 

Note that if r and r' are tilings, we have H T (x) — H T >(x) € Z((l, 1, 1)) and 
h T (x) — h T >(x) e (3), for every x e A . 

Let B = [a|b] be a basis of A with respect to Ao (see Equation (|1|). We say 
that the pair 5(r, B) = (e T (a), e T (b)) composed of the height increment of the 
basic vectors is the fingerprint of r in basis B. We note that for x € Ao and 

y e A, 

h T (x + y) = M x ) + h> T (y), 

and therefore, 

h T (x + [u|v]B(Ai, A 2 ) 4 ) = h T (x) + Aie T (a) + A 2 e T (b), VAi, A 2 e Z. 

In other words, the fingerprint of a tiling determines the height of every point 
in the lattice A. 

3 The Permanent-Determinant method 

In this article we focus the problem of enumerating the A-periodic tilings, for 
a given full-rank sublattice of Ao. Indeed, these tilings correspond to perfect 
matchings in a certain "honeycomb-like" graph (see Figure |8|, obtained as fol- 
lows. We consider two sets of representatives of Ao/A such that the set of 
associated A and y triangles is connected. This way, from A we build (in 
principle not in a unique way) a graph whose vertices are the triangles in a fun- 
damental region: V = A /A x Ao/A. An edge (weighted if desired by L,D,R) 
joins a pair of triangles which can be merged in a lozenge. 
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The Hafnian of a symmetric square matrix of order 2n is defined by: 



Hf (^)= e n a 



mGU-n {ij}6m 



where U n = {{{ii,ji}, ■ ■ ■ , {i n Jn}} I Ui<;<„ {it,ji} = {1, ... ,2n}} is the set of 
the (2n)!/(n!2 n ) (unordered) matchings of the set {1, . . . , 2n}. The number of 
perfect matchings in an undirected graph with an even number of vertices (the 
problem keeps no interest if the number of vertices is odd) equals the Hafnian 
of its adjacency matrix, but the computation of a Hafnian is unfortunately a 
#P-complete problem [16 . The Pfaffian is a related, but tractable, function 
defined over antisymmetric matrices: 



where the parity o~(m) of a matching m = {{i±, ji}, . . . , {i n , in}}) written in 
such a way that ii < ji, for 1 < I < n, is that of the permutation 



It can be shown that, for every antisymmetric matrix A' of even order, we have 
det(-A') = Pi(A') 2 . This result also holds for matrices of odd order, defining 
their Pfaffian as zero. 

Given a directed graph, the Pfaffian of its adjacency matrix A' counts its 
perfect matchings, but some of them affected by a negative sign, leading to an 
(in principle) meaningless sum. Kasteleyn proved (see |7J) that every planar 
graph can be oriented in such a way that Hf(A) = ±Pf(A'), counting therefore 
the perfect matchings. This Pfaffian orientation can be achieved requiring that 
in every face, the number of border edges in clockwise direction is odd. 

However, the graph of Figure [8] is not planar. In general, the graphs we are 
interested in can be embedded in a torus. In [6 , the perfect matching enumer- 
ation problem is solved for the rectangular lattice on a torus, by computing a 
linear combination of four Pfaffians. In a more general fashion, it is stated in 
[8 that the problem can be solved with 4 9 Pfaffians for any graph drawn in a 
surface of genus g. This statement is proved in [5 ; and independently in [14 , 
which contains a general method that uses 2 2 ~ x Pfaffians for a graph embedded 
in a surface of Euler characteristic x, improving therefore the previous state- 
ment for non-orientable surfaces. Tesler method 14J starts drawing the graph in 
a surface represented by a polygon with pasted borders, distinguishing between 
edges contained in the interior of the polygon and those crossing its borders. If 
some of the interior edges form a cycle enclosing the rest, a crossing orientation 
is given to the graph according to the following rule (R4): 

• The set of interior edges is given a Pfaffian orientation, as a planar graph. 



Pf(A') = £ (-!)'(»») J] A\ 
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• Each of the remaining edges is oriented in such a way that any face it 
forms with the interior edges has an odd number of clockwise edges as 
well. 

Before proceeding, let us remark that the graphs we consider are bipartite 
(A may only match y , and conversely). For bipartite graphs, the Permanent- 
Determinant variant, introduced in [12 , simplifies the Hafnian-Pfafhan method. 
Indeed, if G and G are, respectively, an undirected bipartite and a directed 
bipartite graph, let A and A' denote their adjacency matrices, and B and B' 
their bipartite adjacency matrices (whose columns and rows represent black and 
white vertices, respectively). Then, Hf(A) = Per(B) and Pf(A') = det(S'). 

Let A be a full-rank sublattice A of Ao. Firstly, we compute the (unique) 
matrix 

a c 
b 



B = 



e z 2 



such that a > 0, < c < b, and A is generated by the columns of [u|v]B (see 
Equation 0). We choose sets of A and V representatives formed in both cases 
by b rows of a triangles, in the way depicted in Figure [8] More explicitly: 

A : {iv + ju \ < i < b, < j < a}, 

V : {(i + l)v + ju | < i < b, < j < a}. 



We label the edges of the graph with L,D, or R, depending of the orienta- 
tion of the corresponding lozenge. Then, the bipartite adjacency matrix of the 
undirected graph, identifying rows with A and columns with y , presents the 
following decomposition in b x b square blocks of order a: 



M = 



X 

Z X 



Z' 



X 



where Z = D • Id a , 



X = 



R 

L R 



L R 



, and Z' = 





D • Id a _ c ' 


D-Id c 





In the extreme case a = 1, block X equals [R + L]; and when 6=1, the block 
decomposition of M collapses to [X + Z']. Now, in order to define a suitable 
orientation, we start by the edges which do not cross the border of the rectangle, 
orienting by the rule A — > y the D and R edges, and conversely the L edges, 
as is shown in the left side of Figure [9| 

In order to apply Rule (R4) for orienting the rest of the edges, we would need 
that those already oriented are enclosed in a cycle, condition that is not fulfilled 
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Figure 9: Crossing orientation 



in our case. However, we can slightly modify (R4) taking into account those D- 
edges which, joined to some interior edges, form a cycle enclosing an odd number 
of vertices (no R- or L-edges present this pathology under our construction). 
Then, we get a crossing orientation if we orient the border crossing edges as 
follows (see Figure [9]): 

• L-edges: A — » y. 

• D-edges: 

— Id Q _ c : A — > xj, if b is odd; and conversely otherwise. 



Id,. 



V — > A, if b is odd; and conversely otherwise. 



Marking edges crossing the diagonal border with uj± and those crossing the 
horizontal one with ui2, the block structure of the bipartite adjacency matrix M' 
of this directed graph remains as in M, substituting the blocks by: Z — DId a ; 



X 



R 



R 



R 



if a > 1, and X = [R + uiL] otherwise; 





(-l) fc+1 w 2 DId a _ c ' 


(-l)Vw 2 DId c 





Z' = 



We obtain, as a corollary of [14, Theorem 5.2]: 



Theorem 4 Let A be a full-rank sublattice o/Aq. If g{oJ\ , 0J2) = det(M') G 
Z[L, D, R][ui, L02], where M' is defined above, the generating function of the 
K-periodic tilings is: 



Z(L,D,R) - - (.g(l,l) +g(l,-l) +g(-l,l) - g(-l,-l)) . 

It would be interesting to derive a "closed formula" for this generating function, 
or at least for the number of A-periodic tilings Z(l, 1, 1). 
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4 Different tiling types 
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The three constant mappings Tl, Td, and tr are A-periodic tilings for any 
sublattice A C Ao. The identification with piles of cubes gives, for each of these 
tilings, a plane (orthogonal to the Z, X, and Y axes respectively; see Figure 
for the latter two). 

Let B £ Z 2x2 be a 2- rank matrix. We call fundamental triangle to the 
triangle whose vertices are the fingerprints in B of the constant tilings. With 
B = [a|b], we have: 

S(t l ,B) = {-a 1 -2a 2 ,-b 1 -2b 2 ), 
6{t d ,B) = (2a 1 +a 2 ,2b 1 +b 2 ), 
S(t r ,B) = (-a 1 + a 2 ,-b 1 +b 2 ). 

Therefore, the area of the fundamental triangle equals (9/2) det B — (9/2) [Ao : 
A]. Note that the vertices of the triangle are in the lattice (3Z) 2 , possibly shifted 
(see Figure 10). Taking as example the lattice A = [u|v]i?Z 2 , where 



B 



2 -2 
2 4 



the fundamental triangle is defined by the points: 

8{t l ,B) = (-6,-6), 6(t d ,B) = (6,0), S(r R ,B) = (0,6). 









tr 










... v 3 














► » 


* — 


TD 


— 


H 

6 / 


3 

-6 - 


i 1 
















TL 













Figure 10: Types of A-periodic tilings. 



In general, for any two points <5i , <^2 £ Z 2 such that 8± — S 2 € (3Z) 2 , there 
exists a third point 63 = —(Si + 5 2 ) defining with the former the fundamental 
triangle associated to a basis. 



11 



Next result identifies the type of a tiling with a point of the fundamental 
triangle. In the following, when speaking of points in the fundamental triangle, 
we restrict ourselves to those in 5(td, B) + (3Z) 2 . 

Theorem 5 Let A C Ao be a 2-rank lattice and t a A-periodic tiling. Writ- 
ing (L,D,R) = t(r, A) for the type of r, its fingerprint in a base B of A has 
triangular coordinates proportional to (L,D,R): 

S(t, B) = —^—(LS{t l ,B) + DS(r D ,B) + R6(r R , B)). 
[Ao : AJ 

Proof. Let B — [a|b] be the considered basis. We denote by o\ and o 2 the 
orders in the group Ao/A of u and v, respectively, and consider the following 
two representations of A /A: 






u 


2u 




(01 - l)u 


Xl 


Xl + u 


xi + 2u 




xi + (01 - l)u 












x ;i _i 


X ;i _! +U 


x u _! + 2u 




x Zl _i + (01 - l)u 






V 


2v 




(o 2 - l)v 


yi 


yi + v 


yi + 2v 




yi + (02 - l)v 












y« 2 -i 


y« 2 -i + v 


yi a -i + 2v 




y; 2 -i + (o 2 - i)v 



where h = [A : A]/o 4 . For i e {1,2},0 < j < h , let L),D), and i?j be the 
number of elements in the jth row of the ith table whose image by r is L, D, 
and R, respectively. We have: 

L)+D]+R)=o l ,^j. 

There exist constants K\ — (n\,n\,n\) — H t (oiu),k 2 = {k\,k 2 ,k 2 ) = H T (o 2 v) 
such that 

= Lj Rj , Ky == K-y = Dj , 

k\ = k\ = -L), 4 = D 2 j +R). 

Therefore, the numbers Dj , L?, ij + Rj, and D 2 + i?| are independent of index 
j. We consider the matrix 



01 




Ml 


-M2 


2 




-a 2 Zi 


ai?2 



Writing ki = k^ 1 + k? + «f, we have (fci, k 2 ) = (e T (a), e T (b))P, and therefore, 

( e T (a) = aiki/01 + a 2 k 2 /o 2 
\ e T (b) = MiM + b 2 k 2 /o 2 , 
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[A : A]e T (a) = (a 1 / 1 (2Z?) - L) - R)) + a 2 l 2 {D) - 2L) + R])) = 

ai (2D -L-R) + a 2 (D - 2L + R), 
[A : A]e T (b) = (Mi(2Dj - L) - R)) + b 2 l 2 (D] - 2L) + R*)) = 

6i(2L> + L + R) + b 2 (D -2L + R). 

The result follows easily. ■ 
For instance, the tiling shown in Figure [4] is A-periodic where A is defined 



in the example from Figure 10 Its type is (2,2,8); and its fingerprint in basis 
B, (0,3). 



According to Figure 10 there are (at most) seven nonconstant types of A- 
periodic tilings. Four of them involve lozenges in the three different orientations. 
Let us prove that every point 5 in the fundamental triangle represents at least 
one tiling. 

Theorem 6 Let B = [a|b] be a 2-rank integer matrix and A = [u|v]£?Z 2 . The 
set of types {t(r, A) | r G T\} coincides with the intersection of the fundamental 
triangle with S(tjj,B) + (3Z) 2 . 

Proof. As a corollary of Theorem [5j every type lies in that intersection. For 
the converse inclusion, let S be a point in the fundamental triangle. According to 
Proposition [5J the "skeleton" of any tiling with fingerprint S is determined. We 
claim that there is a tiling t with the "skeleton" determined by 5. More formally, 
there is a A-periodic tiling r such that the three-dimensional coordinates of the 
points aiu + a 2 v and &iu + 6 2 v are, respectively, 

A(t, B) = (5i - 2ai — a,2, Si + «i — a 2 , 8\ + a\ + 2a 2 ) , 

- {S 2 - 2b t - b 2 , S 2 + bx - b 2 , 5 2 + bi + 2b 2 ) 
o 

A particular tiling can be formed just by placing a cube under every point of that 
skeleton and filling the position whose coordinates are not bigger component- 
wisely. For instance, applying this process on the skeleton ((—2, 0, 2), (1,-1, 3)) 
we get the tiling depicted in Figure [4] We just need to show that for every 
x, y € A, H T (x) — H T (y) $ (N\{0}) 3 . In other words, that no point in the 
skeleton is "hidden". It is sufficient to check that the normal to the plane 
defined by the two components of A(r, B) lies in the cone {±(x, y, z) £ M 3 | x > 
0,2/ > 0,z > 0}. This condition holds, because of Theorem [5] and the fact that 
the cone is a convex set. I 
The number of points in the border of the fundamental triangle are: 

t d ,t l : gcd(oi + a 2 ,bi + b 2 ) + 1 
t Li t r : gcd(o 2 ,6 2 ) + l 
t r ,t d : gcd(ai,6i) + 1 

By Pick's Theorem, the number of interior points in the fundamental triangle 
equals 

-[A : A] - -(gcd(ai,&i) +gcd(a 2 ,&2) + gcd(ai + a 2l b 1 + b 2 )) + 1 
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and the number of monomials in the generating function Z(L, D, R) is 
^[A : A] + -(gcd(a 1: i»i) + gcd(a 2 , b 2 ) + gcd(ai + a 2) 6i + b 2 )) + 1. 

5 Grouping similar tilings 

The coefficients in the monomials corresponding to borders of the triangle (i.e. 
types involving just one or two lozenge orientations) are easily determined. Let 
us show the case of tiling types with no "L" lozenge. Put d = gcd(ai, bi). There 
are d + 1 points in the edge limited by 6(td,B) and 5(tr, B). These are: 

Si = l -8{T D ,B)+ d ^8(T R ,B), i = 0,...,d. 

According with Theorem[5j a A-periodic tiling r such that 5(r, B) = Si has type 
(i[Ao : A]/d,0, (d — i)[Ao : A]/d). This kind of tilings are constant on the lines 
directed by v and can be enumerated as follows: 

• Compute a triangular form of B: 

'dO' 

• Select % elements x € {0, u, 2u, ...,(d— l)u}, and define r(x) = D for 
them. For the rest, set r(x) = R. 

Therefore, the number of tilings of type (i[A : A]/d, 0, (d - i)[A : A]/d) is (f) 
and there are 2 d tilings whose type is in the considered edge. 

Continuing with the example from Figure [Toj let us enumerate the possible 
tilings in with no "L" lozenge. We need to choose between D and R for r(0) 
and t(u). We get the four tilings depicted in Figure [TT] 

Indeed, the two tilings of type Si = (0,6,6) (with fingerprint (3,3)) are the 
same modulo a shift. In general, we might also be interested in enumerating 
classes of tilings, identifying those which only differ on a shift. This is, defining 
the following equivalence relation on the set of A-periodic tilings T\: 

rSp 3s e A Vx e A : t(x) = p(x + s), 

we need to count elements of T\/S. In the degenerated case of a triangle border, 
the number of classes associated to the point A, is 




the number of necklaces with d beans, i of them coloured. The total number of 

classes of tilings modulo shifts in an edge is - (p(k)2 d ^ k . 

d 

k\d 
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Figure 11: Tilings with only "D" and "R" lozenges. 



In this way, we can define a generating function Z\{L 1 D, R) £ Z[L, D, R] of 
T\/S, analogue to Z(L, D, R), with the same set of monomials indeed. 

Let us consider now an involution which associates pairs of tilings with the 
same type. It is easily derived from Definition |T| that a A-periodic tiling can 
be defined through a bijection f in Ao such that f(x) — x € {0,v — u,v} and 
f (x + y) = t(x) + y, for all x G Ao, y G A. We define then the involution / as 
follows: 

(Jf)(x) = -f- 1 (-x). 

This operation is compatible with the relation £ defined above and it is indeed 
more natural to consider / acting on T^/S. In Section [2j we set the convention 
that the left image in Figure [6] represents a solic cube. Considering the inverse 
convention corresponds to looking at the pile of cubes "from behind". This 
change of viewpoint is encoded by the involution /. Another interpretation 
arises from rotating 180° the plane representation of the tiling. 

We may also find redundant to compute as different tilings related by this 
involution. This allows another simplification in the set T\/S, defining a coarser 
partition. In the degenerated case, the cardinality of the new quotient is the 
number of reversible necklaces with d beads, i of them coloured, but as before, 
the computation of the generating function Zi(L, D, R) seems a more difficult 
problem. 

As we have seen, for types lying in the border of the fundamental triangle 
these functions correspond to well-known combinatorics formulas. However, we 
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are not able to efficiently compute them in the more interesting case of points 
interior to the triangle. 

Let t be a A-periodic tiling containing at least one lozenge in each orienta- 
tions (i.e., no component in its type is zero). This is equivalent to the condition 
of having an inner corner (see Figure [6]) in the tiling. We call flip to a transfor- 
mation of one tiling into another by changing a inner corner into a solic cube, 
or vice versa. It is clear that this operation keeps the type of a tiling. 

For example, a flip in the dotted quoin in the tiling from Figure [4] (indeed, 
that tiling has only a quoin and an inner corner) gives the tiling depicted in 
Figure [T2| which has two quoins and two inner corners. 




Figure 12: 5(t, B) = (0, 3) 



It follows from the proof of Theorem [6] that starting with a A-periodic tiling, 
one may obtain all the tilings with the same type by means of flips. Perhaps 
this consideration is useful in the task on determining Z\ , Z2 ■ On the other 
hand, it is likely that the Hafnian-Pfaffian method, which has proved useful in 
enumerating symmetry classes of tilings (see [10 ), facilitates the evaluation of 
those functions. 
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